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We measured the production of 57 Co, 54 Mn, 68 Ge, 65 Zn, and 60 Co in a sample of Ge enriched in 
isotope 76 due to high-energy neutron interactions. These isotopes, especially 68 Ge, are critical in 
understanding background in Ge detectors used for double-beta decay experiments. They are pro- 
duced by cosmogenic-neutron interactions in the detectors while they reside on the Earth's surface. 
These production rates were measured at neutron energies of a few hundred MeV. We compared 
the measured production to that predicted by cross-section calculations based on CEM03.02. The 
cross section calculations over-predict our measurements by approximately a factor of three depend- 
ing on isotope. We then use the measured cosmic-ray neutron flux, our measurements, and the 
CEM03.02 cross sections to predict the cosmogenic production rate of these isotopes. The uncer- 
tainty in extrapolating the cross section model to higher energies dominates the total uncertainty 
in the cosmogenic production rate. 

PACS numbers: 23.40.-s, 25.40. Fq, 25.40.Sc 



I. INTRODUCTION 

Neutrinoless double-beta decay (Oz//3/3) plays a key role 
in understanding the neutrino's absolute mass scale and 
particle-antiparticle nature 1-6 . If this nuclear decay 
process exists, one would observe a mono- energetic line 
originating from a material containing an isotope sub- 
ject to this decay mode. One such isotope that may un- 
dergo this decay is 76 Ge. Germanium-diode detectors 
fabricated from material enriched in 76 Ge have estab- 
lished the best half-life limits and the most restrictive 
constraints on the effective Majorana mass for the neu- 
trino [7J [5] . One analysis [5j of the data in Ref. [5] claims 
evidence for the decay with a half-life of 2.231q ^ x 10 25 
y. Planned Ge-based 0v(3(3 experiments, MAJORANA [TOl - 
IT2"] and GERDA [T3], will test this claim. Eventually, 
these future experiments target a sensitivity of >10 27 
y or ~1 event /ton- year to explore neutrino mass values 
near that indicated by the atmospheric neutrino oscilla- 
tion results. 

The key to these experiments lies in the ability to re- 
duce intrinsic radioactive background to unprecedented 
levels and to adequately shield the detectors from ex- 
ternal sources of radioactivity. Previous experiments' 
limiting backgrounds have been trace levels of natural 
decay chain isotopes within the detector and shielding 
components. The 7-ray emissions from these isotopes 
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can deposit energy in the Ge detectors producing a con- 
tinuum, which may overwhelm the potential OV/3/3 sig- 
nal peak at 2039 keV. Great progress has been made in 
identifying the location and origin of this contamination, 
and future efforts will substantially reduce this contri- 
bution to the background. The background level goal 
of 1 event /ton-year, however, is an ambitious factor of 
~400 improvement over the currently best achieved back- 
ground level [5]. If the efforts to reduce the natural de- 
cay chain isotopes are successful, previously unimportant 
components of the background must be understood and 
eliminated. The contribution from long-lived isotopes 
produced by cosmic-ray neutrons in Ge detectors fab- 
ricated from enriched Ge was recognized and described 
in [TU US]- In fact, the dominant background that the 
Majorana and GERDA experiments will face, without 
sophisticated analysis cuts, will originate from such iso- 
topes unless mitigation strategies to reduce the activation 
are successful. 

To successfully mitigate the impact of these isotopes 
requires an understanding of their production. Tables [T] 
and |H] summarize the previous production rate estimates. 
The two most critical cosmogenic isotopes, 60 Co and es- 
pecially 68 Ge, show significant variation in the predicted 
rates: factors of ~2 and ~10 respectively. For 60 Co, some 
models predict a higher production rate for enriched Ge 
( enr Ge) than for Ge samples with natural isotopic con- 
tent ( na ^Ge). Reference [16] gives a nice summary of 
previous attempts to calculate the production rates of 
the problematic isotopes, and provides an estimate of its 
own. The authors of that report noted that the calcu- 
lations differ significantly and that measurements would 
be useful to better understand the rates. Some measure- 
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merits do exist, but are either for proton reactions [TT1 - 
[19] , or have a large uncertainty [15] . Since the production 
rates due to neutrons are much larger than for protons, 
its important to have neutron reaction measurements. In 
addition, Barabanov et al. |20] studied the reduction of 
cosmogenic activation as it depends on shielding in order 
to design an optimum transport container. That refer- 
ence also calculated the rates for cosmic-ray proton reac- 
tions and found the production rates for the troublesome 
isotopes to be about a factor of 10 below that for cosmic- 
ray neutrons. The work of Mei et al. [UJ noted that much 
of the large variation in these rates was due to the use 
of different cosmic ray neutron flux estimates, and that 
many of the analyses ( pH El GO] E2 [23] ) used historical 
flux spectra that are less precise than modern measure- 
ments. Although Ref. [21] performed calculations for the 
production of the isotopes of interest to this paper, the 
cross sections were calculated using TALYS [24 . This 
code only predicts cross sections to an energy of 250 MeV, 
whereas other treatments go to higher energies. Refer- 
ence [16] used a modern interpretation of the old neutron 
flux values [25] but not the results of recent measure- 
ments [26] . Furthermore, it used a combination of cross 
section calculations in order to span the energies neces- 
sary for the calculations. Reference [27] presented num- 
bers for a specific shielding geometry, which is not easy 
to translate to a raw production rate with the provided 
data. Hence none of the presently available estimates are 
sufficient to reliably predict the cosmogenic production 
rates and new measurements/estimates are required. 

We exposed a sample of Ge enriched in isotope 76 to 
a wide-band neutron beam that resembles the cosmic- 
ray neutron flux. After exposure we counted the sample 
in a low-background counting system to observe the 7 
rays from the decays of the problematic isotopes. From 
these data we have measured the production rate due to 
fast neutrons in a Ge sample enriched in isotope 76 that 
was taken from material used for Ge detector produc- 
tion. With knowledge of the neutron-beam and cosmic- 
neutron energy spectra, we used these data to provide 
an estimate of the production rate due to exposure of Ge 
to cosmic rays. We used a cross-section calculation that 
spans the energy range of interest and the most recent 
cosmic-ray neutron flux measurements of which we are 
aware. This article describes our determination of values 
for the production rate of these isotopes. 



II. EXPERIMENT 

The sample was exposed to the neutron beam at the 
Los Alamos Neutron Science Center (LANSCE) Weapons 
Neutron Research (WNR) facility from Target 4 Flight 
Path 60 Right (4FP60R) [55] • As the broad-spectrum, 
pulsed neutron beam strikes the Ge target, the outgoing 
7 rays are detected by the GErmanium Array for Neutron 
Induced Excitations (GEANIE) spectrometer [25]. The 
corresponding data from the GEANIE spectrometer will 



TABLE I: A summary of previous estimates of the produc- 
tion of long-lived cosmogenic isotopes in na ^Ge for the iso- 
topes studied in this work. The production rates are given in 
atoms/ (kg d). The data in Ref. [57] was quoted in /xBq/kg 
and we converted to units presented. 



Isotope 


Ref. 
[14] 


Ref. 
[TS] 
(Calc.) 


Ref. 

[TS] 
(Expt.) 


Ref. 

mi 


Ref. 

m 


Ref. 
116] 


Ref. 

\m 


Ref. 
[21] 


57 Co 


0.5 


4.4 


2.9±0.4 


10.2 




9.7 


6.7 


13.5 


54 Mn 




2.7 


3.3L0.8 


9.1 




7.2 




2.7 


68 Ge 


26.5 


29.6 


30±7 


58.4 


82.8 


89 


45.8 


41.3 


65 Zn 


30.0 


34.4 


38±6 


79.0 




77 


29.0 


37.1 


60 Co 


4.8 






6.6 


2.9 


4.8 


2.8 


2.0 



TABLE II: A summary of previous estimates of the produc- 
tion of long-live cosmogenic isotopes in enr Ge for the isotopes 
studied in this work. For these estimates, the abundance val- 
ues are 14% for 74 Ge, 86% for 76 Ge and zero for the other 
naturally occurring isotopes. The production rates are given 
in atoms/ (kg d). 



Isotope 


Ref. 
[14] 


Ref. 
[15] 


Ref. 
|22] 


Ref. 
|20| 


Ref. 

EU 


Ref. 


Ref. 
[211 


This 
Work 


bY Co 


0.1 


1.0 


1.6 




2.3 


2.9 


6.7 


0.7L0.4 


54 Mn 




1.4 


2.3 




5.4 


2.2 


0.87 


2.0L1.0 


68 Ge 


1.2 


1.2 




5.7 


13 


7.6 


7.2 


2.1L0.4 


65 Zn 


6.0 


6.4 


11.0 




24 


10.4 


20.0 


8.9L2.5 


60 Co 


3.5 






3.3 


6.7 


2.4 


1.6 


2.5±1.2 



be used for (n,n'7) analysis that will be presented in a 
separate publication. The GEANIE sample is located a 
distance of 20.34 m from the natural tungsten spallation 
target. 

The target sample was a 11.13-gm, 22-mm diameter 
metal enriched Ge ( enr Ge) powder contained within a 
plastic enclosure. The isotopic abundances within the 
sample were measured by time of flight secondary ion 
mass Spectrometry (ToF SIMS) with the result: 70 Ge 
0.77 ±0.04%, 72 Ge 0.94 ± 0.05%, 73 Ge 0.36 ± 0.03%, 74 Ge 
13.81 ± 0.18%,and 76 Ge 84.12 ± 0.23%. The sample was 
exposed with two separate beam collimations (3/4" and 
1/2" collimators) during 3 irradiation periods. For the 
3/4" -collimator (1/2" -collimator) run a surface area of 
2.85 (1.27) cm 2 was exposed to the beam. The 3/4"- 
collimator exposure was performed between July 16 and 
Juy 23, 2007 (6.99 d elapsed time). The 1/2" -collimator 
exposure was performed between July 27 and August 2, 
2007 (5.95 d) and then between August 8 and August 
14, 2007 (6.19 d). As seen in Fig. [TJ the energy spectrum 
of the third exposure was slightly softer than the other 
2 exposures at the higher energies. The pulsed neutron 
beam has the following timing structure. Macropulses, 
lasting 625 fj,s, occur at a rate of 40 Hz. Micropulses are 
spaced every 1.8 fis, during which the neutron energy 
is determined by the time of flight from the micropulse 
start. An in-beam fission chamber measures the neutron 
flux with 238 U foils. If the reader wishes to convolve the 
neutron spectrum with his/her own cross section model, 
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FIG. 1: The energy spectrum for the neutron beam fluence 
at 4FP60R for the 3 separate exposure periods. The data are 
corrected for live time of the fission chamber. 



FIG. 2: The energy spectrum of 7 rays from the Ge sample 
as measured by a Ge detector (upper curve). Also shown is 
a background spectrum taken with no sample present, (lower 
curve) 



we give a parameterization of the neutron spectrum im- 
pinging upon our sample for convenience. The spectrum 
can be described as: 



III. MEASURED PRODUCTION RATES 



$(E) = (1.325 x 10 



-10\ 



(1) 



The peaks in Fig. [2] were fit to determine the measured 
counts (C) and the results are given in Table III The ef- 



(4.986 In £-3.825 In 2 £+0.9159 In 3 £-0.07402 In 4 E) 

where <& is in units of neutrons/MeV and the energy (E) 
is in MeV. 

The powder was stored for an extended period and 
therefore any radioactivity had decayed to an extremely 
low level (<150 Bq) before counting began. Therefore, 
this sample was well below any action levels and not sub- 
ject to any source-handling requirements. The sample 
was transported to our low-background counting facility 
underground at the Waste Isolation Pilot Plant (WIPP) 
near Carlsbad, NM and counted with a Ge detector. The 
detector was fabricated in 1985 and placed underground 
at WIPP in 1998. It is an n-type semi-coax design with 
a height of 41 mm and a diameter of 51 mm. It is con- 
tained within an «l-mm thick Cu cryostat. The shield 
during these runs consisted of 5 cm of oxygen-free, high- 
conductivity Cu and 10 cm of Pb. 

The sample was counted over a period of 73.86 days 
between February 19 and May 4, 2009 with a total live 
time of 49.02 days. The spectrum is shown in Fig. [2] The 
detector at WIPP has been underground there for over 10 
years and therefore any activitives of isotopes of interest 
to this study, that may have been produced while that 
detector was exposed to cosmic rays as it resided on the 
Earth's surface, have long since decayed away. The lone 
exception to this is a very low level of 60 Co that resides 
in the Cu cryostat of the Ge detector and the inner layer 
of the shield that is also Cu. The background spectrum 
in Fig. [2] shows that this rate is very small compared to 
the sample's 60 Co rate. The data presented in Table III 
includes a subtraction of this background (»T0%) for the 
60 Co count rates. 



ficiency for counting the 7 rays in the geometry used for 
the sample was determined by using 57 ' 60 Co, 54 Mn, 22 Na 
and 137 Cs sources. The geometry of these sources is co- 
incidently very similar to the Ge sample and therefore 
there was no need for simulation to determine the effi- 
ciencies. These sources had calibrated activities known 
to ±1%. Since many of the lines of interest come from 
these isotopes, the measurements were a direct calibra- 
tion of the efficiencies of interest. For 65 Zn and 68 Ga 
lines, a 7-ray efficiency curve was determined from the 
various lines in these sources in order to interpolate be- 
tween the measured line energies. Corrections for sum- 
ming of coincident transitions and annihilation 7 rays 
within the sources were made. The resulting efficiencies 
and their uncertainty resulting from counting statistics 
and the i nter polation and summing corrections are given 
in Table III For the 60 Co calibration, a high event rate 



resulted in a minor peak distortion that the other data 
sets did not suffer. We took an additional uncertainty for 
those peaks to reflect that. We used a random pulser to 
verify the event-rate dependence of the data-acquisition 
system dead time. 

In addition to the 7-ray detection efficiency (e 7 ), an ef- 
ficiency factor must also be included to take into account 
the decay of the isotope since the end of exposure and 
during the counting period. This latter efficiency factor 
(e c ) depends on the half-life of the isotope of interest and 



is also given in Table III 



With these efficiencies, it is straight-forward to cal- 
culate the number of atoms of each isotope that were 
present at the end of exposure to neutrons (August 14, 
2007). The measured number of atoms at the reference 
date of the end of exposure (A*o) is given by: 
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iVn = 



c 



e 7 e c B 



where 



-at:, 



,-AT ! , 



(2) 



(3) 



and T\ nd (Tj tart ) is the number of days since the end 
of exposure that the counting stopped (started) for each 
of the i data runs, A is the decay rate of the isotope in 
question, B is the fraction of decays that emit the 7 ray 
of interest. 

The predicted number of atoms NQ Ted can also be ex- 
pressed in terms of the reaction rates during the 3 expo- 
sure periods: 



JV, 



Pred 



-AT? 



A 



-AT? 



(4) 



where we have corrected for the decay during the expo- 
sure and the decay after the exposure until the reference 
date. In Eqn. [I| T- rrad is the duration time of irradiation 
of the i th exposure, and T de is the time between the 
end of exposure i and the reference date for Nq. ki is 
the production rate (atoms/day) during exposure i and 
is given by: 



A number of systematic effects add to the uncertainty 
in No- These include the uncertainty in the nuclear 
physics parameters, which come from the National Nu- 
clear Data Center. The values for the half-lives and 
branching ratios are known to high precision and are a 
negligible contribution to the total uncertainty. The un- 
certainly in e 7 is described above and included in the 
quoted values. However, each of the sources have an ad- 
ditional uncertainty due to the precision of the known 
activity of 1%. The start and stop times of counting and 
the live time of the counting are known to a small precen- 
tage and are negligible contributions to the uncertainty. 
This is similar for the times associated with the irradia- 
tion. The sample had been stored on the Earth's surface 
for many years prior to exposure to the beam and then 
counting at WIPP. Any isotopes produced by cosmic ray 
neutrons would certainly be at saturation after this ex- 
tended period. The saturation production rates however 
are predicted to be near 1 atom/ (day kg). Therefore the 
total count rate for our 11 g sample due to this contribu- 
tion is several orders of magnitude less than our measured 
count rate and we ignore this systematic effect. 

One can assign an effective energy at which these mea- 
surements were made as that for which the theoretical 
production rate (the product of the CEM cross section 
and the 4FP60R flux) peaks. Table IX provides the ef- 
fective energy for each of the isotopes in question. The 
measured uncertainty in the production rates at these 
energies are given in Table |XI| in the subtotal column. 



hi — 



Fi(E)a(E)dE 



( M S N A 
\{MW)AU 



fi 



(5) 



where n s is the number of exposed sample nuclei per 
unit area in exposure i of duration ti, Fi(E) is the en- 
ergy (E) dependent neutron fluence (number of neutrons 
per MeV) impinging upon the sample during exposure i 
and cr{E) is the energy dependent cross section calculated 
for our sample and its measured isotopic abundances, fi 
is the fluence and cross section energy-dependent inte- 
gral for each exposure. The areal density, n s , is related 
to the total sample mass (M s , 11.13 g) and the sample 
area (A, 3.80 cm 2 ). The beam spot covers most, but 
not all, of the sample geometry and therefore we have ig- 
nored any uncertainty associated with non-uniformities 
in the powdered sample thickness. That is, any such 
non-uniformities will average to the nominal areal den- 
sity. 

]\jPred can wr itt en: 



where is Avagadro's number, and MW is the molec- 
ular weight (75.62 u for our sample). 



IV. CROSS SECTIONS 

If the energy dependence of the neutron fluence over 
the duration of the exposure had been constant and 
identical in shape to the cosmic-ray neutron flux, we 
could have used a simple neutron-flux scaling between 
our measurements and the cosmic-neutron flux to esti- 
mate the cosmogenic production of these isotopes. Since 
the shapes do differ, we used a model for the cross sec- 
tion to adjust for these spectral differences in order to 
make a prediction regarding the cosmic-ray production 
rates. For this purpose, we use cross sections for iso- 
tope production calculated using Cascade-Exciton Model 
(CEM03.02 (HH33) as it usually has a better predictive 
power in comparison with other similar available models 
(see, e.g. [33]). The CEM formalism permits the calcu- 
lation of cross sections to the high energies necessary for 
this work. We assumed that the energy dependence of 
the cross section is nearly correct but that there may 
be an overall normalization uncertainty. This procedure 
then corrects for any normalization uncertainty. Below 
we describe how we estimate the uncertainty associated 
with the energy dependence assumption. 

To determine the yield of an isotope, one must consider 
all feeder isotopes that may decay to the isotope of inter- 
est. Therefore to determine the cumulative cross section 
for 57 Co for example, one must sum the cross sections 
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TABLE III: A summary of No for the long-lived isotopes in the Ge sample. The uncertainty quoted for JVo is statistical and 
arises only from the counting statistics of the measured peaks. The quoted uncertainty for N Pred arises from the uncertainty 
in the 4FP60R neutron fluence. Systematic uncertainties are summarized in Table |Xl] and discussed in the text. 



Isotope 


n/2 (days) 


7-ray Energy 


C 


e 7 


e c 


B 


iVo 


N gred 


bv Co 


271.8 


122.1 keV 


2916±84 


0.1663(7) 


0.0280 


0.856 


7.31 ± 0.21 x 10 b 


2.97 ±0.06 x 10 b 


57 Co 


271.8 


136.5 keV 


386±63 


0.163(2) 


0.0280 


0.107 


7.89 ± 1.3 x 10 5 


2.97 ±0.06 x 10 6 


54 Mn 


312.1 


834.9 keV 


1084±43 


0.0302(3) 


0.0297 


1.000 


1.21 ± 0.05 x 10 6 


1.79 ±0.05 x 10 6 


68 Ge 


270.8 


1077.4 keV 


198±18 


0.0207(5) 


0.0280 


0.032 


1.06 ± 0.10 x 10 7 


2.92 ±0.02 x 10 7 


65 Zn 


244.25 


1115.5 keV 


8541±95 


0.0232(5) 


0.0262 


0.506 


2.77 ±0.03 x 10 7 


5.03 ±0.06 x 10 7 


60 Co 


1923.6 


1173.2 keV 


1342±42 


0.0200(6) 


0.0146 


0.999 


4.61 ± 0.14 x 10 6 


7.50 ±0.01 x 10 6 


60 Co 


1923.6 


1332.5 keV 


1176±39 


0.0167(3) 


0.0146 


1.000 


4.83 ± 0.16 x 10 6 


7.50 ±0.01 x 10 6 



TABLE IV: The calculated cumulative cross sections for pro- 
duction of selected cosmogenic isotopes in 70 Ge. The cross 
sections are cumulative because they are sums over all iso- 
topes that feed the isotope of interest. 



Energy 
(MeV) 


Isotope Production 


Cross Section (mb) 


b 'Co 


°"Mn 


bs Ge 


00 Zn 


bu Co 


10 

















20 

















30 








295.6 


16.63 





40 








248.0 


97.52 





50 








153.1 


51.05 


0.0022 


60 








119.6 


36.86 


0.034 


70 








101.6 


55.99 


0.53 


80 


0.009 


0.001 


90.53 


86.66 


0.80 


90 


0.019 


0.0019 


82.67 


91.14 


0.70 


100 


0.16 


0.0037 


75.15 


88.94 


0.80 


200 


11.49 


3.06 


50.42 


54.34 


4.78 


300 


16.52 


7.72 


43.50 


47.17 


4.34 


400 


18.97 


10.59 


38.13 


41.63 


4.38 


500 


20.34 


12.86 


33.34 


37.14 


4.36 


600 


20.28 


14.31 


29.17 


33.39 


4.21 


700 


19.09 


14.49 


26.46 


30.89 


3.81 


800 


18.59 


14.15 


24.75 


28.58 


3.63 


900 


17.28 


13.59 


23.36 


27.026 


3.34 


1000 


16.27 


13.08 


22.07 


25.99 


3.11 


2000 


10.22 


8.12 


16.62 


20.61 


2.17 



TABLE V: The calculated cumulative cross sections for pro- 
duction of selected cosmogenic isotopes in 72 Ge. The cross 
sections are cumulative because they are sums over all iso- 
topes that feed the isotope of interest. 



Energy 
(MeV) 


Isotope Production Cross Section (mb) 


b 'Co 


M Mn 


b8 Ge 


bb Zn 


bu Co 



















































































15.82 


0.15 














45.86 


13.43 














42.57 


20.87 


0.0010 











34.56 


18.01 


0.014 











32.2 


19.43 


0.20 











30.69 


30.27 


0.61 


4.17 


4.32 


0.99 


22.9 


42.50 


5.26 


9.45 


10.00 


5.03 


19.88 


38.31 


5.78 


12.06 


12.86 


8.55 


16.91 


34.84 


6.22 


14.25 


15.12 


11.47 


14.77 


31.37 


6.56 


15.13 


16.05 


13.91 


13.46 


28.30 


6.45 


14.97 


15.95 


14.48 


12.45 


25.54 


6.09 


14.29 


15.22 


14.62 


11.54 


23.63 


5.85 


13.46 


14.37 


14.54 


11.13 


21.75 


5.52 


12.62 


13.45 


14.07 


10.65 


20.59 


5.14 


7.72 


8.24 


8.79 


7.96 


14.69 


3.27 



for 57 Co, 57 Ni, 57 Cu and 57 Zn. Similarly for 68 Ge, one 
must also consider 68 As and for 65 Zn one must consider 
65 Zn, 65 Ga, 65 Ge and 65 As. For the isotopes 54 Mn and 
60 Co one only need consider the primary isotopes. Note, 
that for the case of 60 Co, this is only approximately true. 
Because of the long half-life of 60 Fe, the contributions of 
60 Fe, 60 Mn and 60 Cr are negligible. Using the CEM3.02 



code, Tables IV through VIII provide the values of these 
cumulative cross sections summed over the feeders for the 
5 isotopes that comprise Ge. 

Table |IX| gives the critical neutron energy range for 
the majority of the isotope production rate. This low 
(high) value of this energy range is defined as the energy 
for which 10% (90)% of the cumulative production rate 
is reached. The ranges are given for both the 4FP60R 
beam spectrum and the cosmic ray flux. The 4FP60R 
range also is calculated for the isotopic abundance of our 
sample, whereas the cosmic ray flux column is calculated 



for our standard 86%-14% isotopic enrichment. Although 
the energy at which the peak of the production rate for 
the two spectra are similar, the cosmic ray production 
rate is still significant at higher energies. Figure [3] shows 
the production rate as a function of energy for the cos- 
mic ray flux. Since there is a spectral difference between 
the 4FP60R flux and the cosmic ray flux above 200 MeV, 
we must use theory along with our measurements to es- 
timate the cosmic-ray production. Therefore, we must 
consider the possibility that an energy dependence in the 
theoretical cross section used for this extrapolation will 
introduce a systematic uncertainty. The input physics to 
the CEM03.02 cross sections does not change over the 
energy range of interest and no known physics is omit- 
ted. Therefore one does not expect the variation in the 
cross section uncertainty with energy to be large. Refer- 
ences [5DH51] indicate that the predictions of the model 
are uniform across that range. Reference |33) presents 
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TABLE VI: The calculated cumulative cross sections for pro- 
duction of selected cosmogenic isotopes in 73 Ge. The cross 
sections are cumulative because they are sums over all iso- 
topes that feed the isotope of interest. 



Energy 
v ivie v ) 


Isotope Production Cross Section (mb) 


uo 


Mil 


i^e 


ZjTL 


Oo 


10 

















20 

















30 

















40 

















50 

















60 








4.79 


0.08 





70 








22.14 


6.67 





80 








23.51 


13.28 


0.002 


90 








19.93 


11.61 


0.017 


100 








18.37 


12.60 


0.13 


200 


3.63 


0.52 


15.22 


35.05 


4.94 


300 


7.45 


3.87 


14.00 


33.30 


6.14 


400 


10.40 


7.34 


11.98 


30.62 


6.88 


500 


12.84 


10.60 


10.39 


27.92 


7.32 


600 


14.07 


13.07 


9.20 


25.63 


7.37 


700 


13.83 


14.45 


8.41 


23.38 


7.17 


800 


13.67 


14.57 


7.84 


21.48 


6.91 


900 


13.06 


14.49 


7.25 


19.88 


6.51 


1000 


12.51 


13.81 


6.92 


18.58 


6.19 


2000 


7.73 


8.64 


5.53 


12.96 


3.94 



TABLE VII: The calculated cumulative cross sections for pro- 
duction of selected cosmogenic isotopes in 74 Ge. The cross 
sections are cumulative because they are sums over all iso- 
topes that feed the isotope of interest. 



Energy 
(MeV) 


Isotope Production Cross Section (mb) 


bY Co 


M Mn 


m Ge 


0S Zn 


""Co 


10 


0.00 


0.00 


0.00 


0.00 


0.00 


20 


0.00 


0.00 


0.00 


0.00 


0.00 


30 


0.00 


0.00 


0.00 


0.00 


0.00 


40 


0.00 


0.00 


0.00 


0.00 


0.00 


50 


0.00 


0.00 


0.00 


0.00 


0.00 


60 


0.00 


0.00 


0.00 


0.00 


0.00 


70 


0.00 


0.00 


0.57 


0.00 


0.00 


80 


0.00 


0.00 


7.78 


1.61 


0.00 


90 


0.00 


0.00 


13.04 


6.79 


0.00 


100 


0.00 


0.00 


12.79 


8.06 


0.01 


200 


1.03 


0.23 


10.76 


28.71 


4.07 


300 


5.37 


2.73 


10.43 


28.24 


6.04 


400 


8.16 


6.14 


8.74 


27.13 


7.28 


500 


10.85 


9.52 


7.69 


25.33 


8.09 


600 


12.16 


12.18 


6.85 


23.42 


8.08 


700 


12.77 


13.66 


6.17 


21.17 


8.19 


800 


12.40 


14.34 


5.44 


19.32 


7.64 


900 


11.96 


14.44 


5.11 


17.64 


7.32 


1000 


11.36 


13.93 


4.78 


16.47 


6.88 


2000 


7.03 


8.76 


3.79 


11.31 


4.46 
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FIG. 3: The CEM cross sections folded with the cosmic-ray 
neutron energy spectrum as a function of energy for the iso- 
topes of interest. 



a quantitative analysis of the comparison between the 
model and previous experiments using a large number of 
data sets. The result is that the uncertainty is fairly en- 
ergy independent with a 2-10% change depending on the 
data set used for comparison. 

Unfortunately there is a lack of cross section data for 
large AA transitions induced by fast neutrons. Such data 
would directly validate the theoretical model for inter- 
preting our production rate measurements in terms of 
the cosmic-ray neutron flux. Ref. [33], however, does 
provide data on proton-induced transitions on an 56 Fe 
target with AA values similar to those of interest in this 
work. (See Figs. 3-5 in that work.) For each isotope con- 



sidered here, we use the data of Ref [33J to estimate the 
uncertainty that may arise from the energy dependence 
of the cross section. 

The critical energy range for the cosmic-ray produc- 
tion rate for the most important /3/3-background isotope, 
68 Ge, is between 110 and 530 MeV. (The low (high) en- 
ergy of this range was defined by integrating over energy 
up to the 10% (90%) value of the total production rate.) 
Ref. [33] studied a few AA=8 transitions in 56 Fe (similar 
to the production of 68 Ge from 76 Ge) that can be used 
to estimate this uncertainty for the production of 68 Ge. 
The AA=8 transitions in 56 Fe producing nuclei 48 Cr, 48 V 
and 48 Sc are then used for this comparison. The ratio of 
the measured and calculated cross sections for these pro- 
cesses varies by an average of ~15% between 300 and 500 
MeV. We estimate the uncertainty in the extrapolation 
to higher energies as being equal to this 15% variation 
scaled to the fraction of model-predicted cosmic-ray pro- 
duction above 200 MeV. The result is given in Table [IX] 

This uncertainty contribution is different for each iso- 
tope due to the different critical neutron energy ranges, 
the different fraction of cosmic-ray production rate over 
200 MeV, and the different results of the comparison of 
theory to data in Ref. [33]. For the AA=11 production 
of 65 Zn, we used the data and theory of the isotope pro- 
duction of 44 Sc and 44 K between 300 and 750 MeV with 
the result that the ratio varies «35%. For the AA=16 
production of 60 Co, we used the data and theory of the 
isotope production of 41 Ar and 39 CI between 300 and 750 
MeV with the result that the ratio varies ~50%. For the 
AA=19 production of 57 Co, we used the data and the- 
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TABLE VIII: The calculated cumulative cross sections for 
production of selected cosmogenic isotopes in 76 Ge. The cross 
sections are cumulative because they are sums over all iso- 
topes that feed the isotope of interest. 



Energy 
(MeV) 


Isotope Production Cross Section (mb) 


bY Co 


M Mn 


M Ge 


BS Zn 


uu Co 


10 


0.00 


0.00 


0.00 


0.00 


0.00 


20 


0.00 


0.00 


0.00 


0.00 


0.00 


30 


0.00 


0.00 


0.00 


0.00 


0.00 


40 


0.00 


0.00 


0.00 


0.00 


0.00 


50 


0.00 


0.00 


0.00 


0.00 


0.00 


60 


0.00 


0.00 


0.00 


0.00 


0.00 


70 


0.00 


0.00 


0.00 


0.00 


0.00 


80 


0.00 


0.00 


0.00 


0.00 


0.00 


90 


0.00 


0.00 


0.04 


0.00 


0.00 


100 


0.00 


0.00 


1.00 


0.19 


0.00 


200 


0.12 


0.11 


5.52 


17.18 


1.95 


300 


2.43 


1.24 


5.92 


19.91 


5.07 


400 


4.92 


3.94 


5.27 


20.67 


7.15 


500 


7.39 


7.15 


4.69 


20.14 


8.52 


600 


9.16 


10.34 


4.11 


18.97 


9.12 


700 


9.84 


12.25 


3.73 


17.41 


9.40 


800 


10.21 


13.63 


3.36 


16.11 


9.16 


900 


10.02 


13.86 


3.08 


14.84 


8.66 


1000 


9.62 


14.00 


2.85 


13.60 


8.36 


2000 


5.82 


8.92 


1.82 


8.67 


5.43 



V. CONVERTING THE MEASURED 
PRODUCTION RATES TO COSMOGENIC 
PRODUCTION RATES 

To estimate the production rate due to cosmic-ray neu- 
trons, one must know the energy dependent flux of the 
neutrons and the cross sections. 

Ziegler carried out a comprehensive study on the 
cosmic-ray neutron flux |25j and pointed out that some 
of the data from early measurements is incorrect or of 
marginal quality. Mei et al. [22 recognized that previous 
estimates of the cosmogenic production rates used vari- 
ous outdated estimates of the cosmic neutron flux. Im- 
proved recent measurements by Gordon et al. |26| show 
that the flux density spectrum at sea level can be param- 
eterized as 

4>{E) = 1.006 X In 2 £+2.1451 In £ 

+ 1.011 X l -3 e -0.41061n^0.6671n£ (?) 

where E is neutron kinetic energy in MeV and <f> is given 
in units of cm _2 s _1 MeV _1 . This parameterization func- 
tion agrees with the data within ~2%. Note that param- 
eterization used by Ref. [16] based on that from Ref. [25] 
differs from that of Ref. [55] . Both curves are shown in 
Fig-! 



ory of the isotope production of 38 C1 and 38 S between 
300 and 1000 MeV with the result that the ratio varies 
«50%. For the AA=22 production of 54 Mn, we used the 
data and theory of the isotope production of 38 CI and 
28 Mg between 300 and 1000 MeV with the result that 
the ratio varies «50%. All these results are summarized 



in Table IX along with the deduced uncertainty in the 
estimated cosmic-ray production rate. The uncertainty 
associated with this scaling dominates the total uncer- 
tainty for the cosmic-ray production rate estimates (see 
Table [Xlj). Because the production rate of 68 Ge is concen- 
trated at lower energies, the theoretical extrapolation for 
this important isotope is more reliable than the others. 



A reader may prefer to use a different cross section 
model to estimate the production due to cosmic rays. If 
so, Eqns. [T] and [7] can be used with a different model 
to interpret our measurements of the production rate at 
WNR in terms of the cosmic ray flux. Furthermore, the 
uncertainty assigned to the choice of the cross-section 
model's energy dependence is not based on neutron pro- 
jectile data, due to a lack of such data. Therefore, one 
must recognize the caveat this introduces in the extrap- 
olation of the measured production rate values to pre- 
dicted cosmic-ray production rates. 




en 
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FIG. 4: The measured neutron flux parameterization func- 
tions at sea level |25U26| and a normalized plot of the 4FP60R 
neutron fluence showing that the spectral shape in the critical 
100-300 MeV range is similar to the cosmic ray spectrum. 

Gordon's measurements show that the shape of the 
outdoor ground-level neutron spectrum does not depend 
significantly on altitude, cutoff or solar modulation. Fac- 
tors depending on atmospheric depth, geomagnetic cut- 
off, rigidity, and geomagnetic location are given by Gor- 
don for correcting the flux for these effects. His model 
is slightly different than Ziegler or the JEDEC standard. 
As a comparison, the cosmic production rates resulting 
from Gordon's flux are 13-33% higher than that using the 
flux of Ziegler for the isotopes of interest in this work. 



TABLE IX: A summary of the estimates of the uncertainty in the cosmic-ray production rate due to the uncertainty in the 
energy dependence of the cross section extrapolation. The second column gives the neutron energy range over which the CEM 
code predicts the majority of the production rate for the 4FP60R beam spectrum and our isotopic sample. The third column 
gives the neutron energy range for which the CEM code predicts the majority of the production rate for a cosmic ray beam 
spectrum and a 86%-14% enriched sample. The fourth column gives the fraction of the production rate arising from neutrons 
with energy less than 200 MeV. The fifth column is the estimated uncertainty in the energy dependence of the cross section 
as described in the text. The last column provides the deduced uncertainty contribution to the cosmic-ray production rate 
resulting from scaling our measurements. 



Isotope 


Critical E n for 
Beam Prod. Rate 


Critical E n 
for Cosmic Prod. 


Fraction Cosmic Prod, 
for E n < 200 MeV 


Est. Uncertainty in Energy 
Dependence of Cross Section 


Deduced Uncertainty in 
Production Rate (%) 


bv Co 


170 - 420 MeV 


280 - 1040 MeV 


1.9% 


50% 


49% 


54 Mn 


200 - 490 MeV 


330 - 1190 MeV 


1.0% 


50% 


50% 


68 Ge 


35 - 250 MeV 


110 - 530 MeV 


44.2% 


15% 


8.4% 


65 Zn 


110 - 300 MeV 


140 - 640 MeV 


29.2% 


35% 


24.8% 


60 Co 


140 - 350 MeV 


190 - 890 MeV 


12.4% 


50% 


43.8% 



TABLE X: A summary of the parameters that enter into the 
calculation of our estimate for the cosmic-ray production rate 
of the various isotopes (K scaled)- Th e uncertainties quoted 

The value for 57 Co is 
has a much smaller 



XI 



for K sca i et i are explained in Table 
taken as that of the 122-kev line as it 
total uncertainty. The 60 Co value is the average of the two 
individual measurements and, since the uncertainty is dom- 
inated by systematic effects, we quote the larger of the two 
associated total uncertainty values given in Table |XI| These 
production rate values are for a nominal 86% 76 Ge 14% 74 Ge 
isotopic mixture. 



Isotope 
Isotope 


Ratio 
f N ° ) 


K Gordon 

(atoms/ (kg d) 


K scaled 

(atoms/ (kg d) 


bv Co 


0.25 ±0.01 


2.93 


0.72 ±0.37 


54 Mn 


0.67 ±0.03 


2.91 


1.96 ± 1.01 


68 Ge 


0.36 ±0.03 


5.83 


2.12 ±0.39 


65 Zn 


0.55 ±0.01 


16.24 


8.94 ±2.53 


60 Co 


0.63 ±0.03 


4.06 


2.55 ± 1.20 



systematic uncertainty of 13.5%. 

We estimated the rate of cosmogenic production in 
enr Ge with isotopic abundance limited to isotopes 74 
(14%) and 76 (86%). However, a small admixture of iso- 
tope 70 can result in significantly more 68 Ge due to the 
much lower threshold for the neutron reaction and the 
higher cross section. Here we present a rough formula 
to estimate the 68 Ge production rate (Pes) a s a function 
of the amount of 70 Ge present in the sample in per cent 
(X). 



P 6S = (2.12 ± 0.75X)&toms/(kg d) 



(8) 



This formula should be accurate to approximately 
20%. The other critical isotope for f3/3 research is 60 Co 
and it has a small production dependence on the fraction 
of 70 Ge. 



The ratio of N to values given in Table 



III 



be used to scale a predicted production rate (K Gordon) 
based on the cross sections of Section IIVI and the Gor- 
don neutron flux parameterization to provide an estimate 
(K sca i e d) of the cosmogenic production rate indicated by 
our measurements. The numbers used in this arithmetic 
are summarized in Table [X] with the results also quoted 
in Table [TT| for comparison to previous predictions. 

The estimate of the cosmogenic production from these 
measured results have some additional systematic uncer- 
tainties. These include the precision to which the cos- 
mic ray neutron flux is known (10-15% and so we split 
the difference and use 12.5% [26J. The total cosmogenic 
rate includes contributions from subdominant proton and 
pion interactions. These only contribute approximately 
10% [5U] to the total production rate. These charged 
particles are much less penetrating than neutrons and 
therefore their impact on any given sample is very ge- 
ometry dependent. Hence we assume a 50% uncertainly 
on this correction for an uncertainty of 5%. Again these 
uncertainties are uncorrelated and result in a estimated 



VI. DISCUSSION AND CONCLUSION 

We measured the production of 57 Co, 54 Mn, 68 Ge, 
65 Zn, and 60 Co in a sample of Ge enriched in isotope 76 
due to high-energy neutron interactions within a neutron 
beam with a spectrum similar to that of the cosmic-ray 
neutron flux at the Earth's surface. The uncertainty on 
the beam-produced production rate at those energies is 



below 10%, depending on isotope, as given in Table XI 



The results, presented in Table |III[ were compared to 
cross sections calculated with CEM03.02. The measure- 
ments are smaller by about a factor of 2-4 than these 
calculations depending on isotope. We scaled these mea- 
surements using CEM03.02 cross sections to estimate the 
cosmic-ray production rate of these troublesome isotopes 
and present the results in Table |Xj Our estimated to- 
tal uncertainty for the production rate of the critical 
68 Ge isotope is approximately 20% providing much bet- 
ter guidance to double-beta decay experimenters in their 
efforts to understand background due to this isotope. 
The uncertainties on the deduced cosmogenic production 
is dominated by the uncertainty in the extrapolation of 
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TABLE XI: A summary of the uncertainties (in %) that contribute to the total uncertainty of cosmic production rate. The 
column labeled SubTotal refers to the quadrature sum of all non-cosmic-ray related contributions and provides the uncertainty 
for the ratio in Table IXl 



T^ntnnp 

i o w U w l_t 


r!miTitincr 


Piffi PI PT1 PV 




P rpH i pt pH 

J. i. l^k-llV^ ICQ 


Flux 


SubTotal 


C^pimrnp 


l\fpiif vnvi 

1 l C LI Ij L l_JLL 


Prot on 


Total 


(Line Energy) 


Statistics 




Activity 


4FP60R 


Chamber 




Neutron 


Specral 


Correction 












Fluence 


Live Time 




Flux 


Difference 






bY Co(122) 


2.9 


0.4 


1.0 


0.9 


0.3 


3.2 


12.5 


49 


5.0 


51 


57 Co(136) 


16.3 


1.2 


1.0 


0.9 


0.3 


16.4 


12.5 


49 


5.0 


53 


54 Mn 


4.0 


1.0 


1.0 


1.2 


0.3 


4.4 


12.5 


50 


5.0 


51 


68 Ge 


9.1 


2.4 


1.0 


0.5 


0.3 


9.5 


12.5 


8.4 


5.0 


18 


65 Zn 


1.1 


2.2 


1.0 


0.5 


0.3 


2.7 


12.5 


24.8 


5.0 


28 


60 Co(1173) 


3.0 


3.0 


1.0 


0.6 


0.3 


4.5 


12.5 


44 


5.0 


46 


60 Co(1332) 


3.2 


1.8 


1.0 


0.6 


0.3 


3.9 


12.5 


44 


5.0 


46 



the cross sections to higher neutron energies. This large 
uncertainty arises due to a lack of large AA, neutron- 
induced reactions at higher incident energies with which 
to validate the theory used for the extrapolation. 
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